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Metal-insulator transitions in pyrochlore iridates (A2lr207) are believed to occur due to subtle interplay of 
spin-orbit coupling, geometric frustration, and electron interactions. In particular, the nature of magnetic order- 
ing of iridium ions in the insulating phase is crucial for understanding of several exotic phases recently proposed 
for these materials. We study the spectrum of magnetic excitations in the intermediate-coupling regime for the 
so-called all-in/all-out magnetic state in pyrochlore iridates with non-magnetic A-site ions (A=Eu,Y), which 
is found to be preferred in previous theoretical studies. We find that the effect of charge fluctuations on the 
spin-waves in this regime leads to strong departure from the lowest-order spin-wave calculations based on mod- 
els obtained in strong-coupling calculations. We discuss the characteristic features of the magnetic excitation 
spectrum that can lead to conclusive identification of the magnetic order in future resonant inelastic x-ray (or 
neutron) scattering experiments. Knowledge of the nature of magnetic order and its low-energy features may 
also provide useful information on the accompanying metal-insulator transition. 



I. INTRODUCTION 

Pyrochlore iridates (A2lr207) have recently attracted much 
attention as prominent examples of 5d transition metal ox- 
ides where interplay of spin-orbit (SO) coupling and elec- 
tron interactions can lead to a number of competing ex- 
otic phases.'"-' Interestingly, most of these materials are 
found to either exhibit a finite temperature metal-insulator 
(MI) transition'"'"'^"'^-"^*', or naturally lie close to a zero 
temperature MI quantum phase transition (both pressure- 
driven^"'-^, and/or chemical -pressure-driven via variation of 
A-site ions'"'-^). It is now beUeved that the nature of the 
MI transitions in these systems are related to the magnetic 
order of iridium (Ir) ions in the low-temperature insulating 
phase. '"•-''"-^ Further, the details of such magnetic ordering 
pattern are known to be crucial for some of the proposed 
novel phases like the Weyl semimetal.^^-^ Therefore, the de- 
termination of Ir magnetic configuration is important, both to 
gauge the relevance of the proposed novel phases and to shed 
light on the nature of the MI transition in pyrochlore iridates. 
However, conclusive experimental evidence for the nature of 
such magnetic order in pyrochlore iridates with non-magnetic 
A-site ions (such as EuiIriO? and Y2lr207)"''2'''''22,23,25,28 
presently lacking. 

Several important clues regarding the magnetic order have 
been revealed by recent muon spin resonance/relaxation 
and magnetization measurements on both Eu2lr207" -^ and 
¥2^207''^: in the low temperature insulating phase, these 
measurements suggest that localized Ir moments exhibit long- 
range magnetic order that may not break the pyrochlore lattice 
symmetry. These results are consistent with the claim that in 
the ground state, the Ir moments order in the non-collinear, 
all-in/all-out (AIAO) fashion (see Fig. 1), as was previously 
found in calculations for both the strong-^ and intermediate 
electron correlation regimes. ^^'^ However, the above findings 
cannot conclusively prove that Eu2lr207/Y2lr207 orders in the 
AIAO fashion; a study of the low-energy magnetic excitations 
is required to identify the signatures unique to the AIAO state. 




FIG. 1. (Color online) The top left tetrahedron indicates the sublat- 
tice indices used in our work. The red vector in the top right tetrahe- 
dron shows the nearest-neighbor d32. The gold vectors pi4 and qi4 
are perpendicular to the C2 axis and span the plane in which next- 
nearest-neighbor V((i4» must lie. The bottom right tetrahedron with 
the green vectors shows the all-in/all-out configuration. 



Indeed, recent RJXS and inelastic neutron scattering experi- 
ments have revealed the magnetic excitation spectra of other 
iridates such as Sr2lr04^^°, Sr3lr207^', and Na2lr03^^-. 

In this paper, we compute the magnetic excitation spec- 
trum in both the intermediate and strong electron correlation 
regimes. Starting from an effective Hubbard model in the 
/efF =1/2 basis to describe Ir electrons, we study the evolu- 
tion of the magnetic excitation spectrum in the intermediate- 
regime by computing the transverse magnetic dynamic struc- 
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ture factor within the random phase approximation (RPA). 
The robust symmetry -protected degeneracies of the spin-wave 
spectrum, the occurrence of Landau damping, as well as the 
characteristic dispersion along high symmetry directions are 
demonstrated to be the defining signatures of the AIAO state 
in the intermediate- 1/ regime. In the strong correlation limit, 
we derive an effective spin-model and calculate the corre- 
sponding spin- wave excitations. Through naive fitting of the 
RPA dynamic structure factor with the strong-coupling spin- 
wave results, we find that the magnetic excitation spectrum at 
intermediate- [/ shows strong departures from the lowest-order 
spin-wave calculations. Current experiments suggest that the 
intermediate-coupling regime may be more suitable for de- 
scribing Eu2lr207, where the charge gap estimated from the 
resistivity measurements is found to be small (~ 10 meV)." 

The rest of this paper is organized as follows. We be- 
gin with a brief description of the Hubbard model relevant 
to the pyrochlore iridates in Section II. We introduce a new 
parametrization for the tight-binding parameters, which eluci- 
dates the structure of the mean-field phase diagram presented 
in Section III. Subsequently, we discuss the results obtained 
both in the intermediate- [/ (Section IV A) and large-[/ (Sec- 
tion IV B) regimes. We also discuss the unique signatures of 
the AIAO phase and compare the results obtained in the two 
different regimes. Implications of our results are summarized 
in Section V. Further details regarding various calculations are 
given in the appendices. 



II. MICROSCOPIC HAMILTONIAN 

In the pyrochlore iridates, Ir"*^ ions are located at the cen- 
ter of corner-sharing oxygen octahedra. This leads to crystal- 
field splitting of the 5d orbitals into upper eg and lower t2y 
orbitals. The six-fold degenerate t2g orbitals (including spin 
degeneracy) is then split into an upper 7eff =1/2 doublet and 
lower Jgff = 3/2 quadruplet by the atomic SO coupling with 
energy separation of 3A/2 (A ^ 500 meV for Ir). Therefore, 
in the atomic Umit, the five valence electrons will fully fill 
the Jgff -3/2 states, half-fill the 7etf =1/2 states, and leave 
the eg orbitals unoccupied. This atomic picture suggests that 
the low-energy physics can be captured by considering only 
the half-filled /eft = 1/2 states.** ** Therefore it is useful to start 
from the most general on-site Hubbard model in the 7eff =1/2 
basis allowed by symmetry: 



(1) 



where i,j are lattice site indices, c,- - (c,|, c,|) are the elec- 
tron annihilation operators, t and i are the z-components 
of the psuedospin operator defined in the global basis, and 
Ilia- = cj^c/o- is the electron number operator at site / of psue- 
dospin cr. In general the hopping matrix hij is complex and 
can be constrained by considering time-reversal invariance 
and various space group symmetries (Moriya rules). 

Time-reversal invariance restricts the hopping matrix to the 
form: 



where I is the 2x2 identity matrix, cr{- a^, (Ty, cr.) are the 
Pauli matrices (in the pseudospin space), and tjj and v,j = 
(v^*^., v'-fj, vjP are real hopping amplitudes. Hermiticity of the 
Hamiltonian impUes 



and 



(3) 



tij and \jj transform as a scalar and as a psuedovector respec- 
tively under the space group symmetries of the lattice.-''''''' 
In particular these symmetries constrain the nearest-neighbor 
(NN) V(,;)S to be perpendicular to the mirror plane contain- 
ing / and j. The two possible directions that v<,j) can point 
are classified as "direct" or "indirect".-'' We shall denote the 
unit vector of the "direct" case as d,j (see Fig. 1). For next- 
nearest-neighbors (NNN) {{ij}}, the two-fold rotation axis that 
exchanges / and j restricts to lie in the plane normal to 
this axis.''"* We can parametrize this plane by the two orthonor- 
mal vectors 



P/;= V6/4(-R,v + Do) 
q,V= V3/2(r,-, + D,v), 



(4) 
(5) 



where R,y = x r^j and s da x d^j, site k being the 
common NN of / and j (see Fig. 1). 

The above constraints allow us to parametrize hij by two 
real NN (unprimed) hopping amplitudes and three NNN 
(primed) hopping amplitudes: 



h(ij){ti,t2) - t\l + itodij ■ cr. 



(6) 



For the rest of the paper, instead of using (f i , to), (f j , t'^, t'^ as 
our hopping parameters, we find it more convenient to work 
with (t,0),(t',6',(p') which makes certain symmetries of the 
phase diagram readily accessible. These are defined as 
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sin(6',/2 - 0) 



sm(e,/2 - 9') sin(0') 
cos(0') 



(7) 



(8) 



hij - tiji + ■ cr. 



(2) 



where 9, - 2arctan(V2) ^ 109.47° is the tetrahedral an- 
gle. This parametrization together with the definition of p,j 
and q,y naturally makes the following property of the model 
manifest: if we perform the basis transformation Cina — > 
e^'™"'°'"i'Ci„j}, where h„ is the unit vector pointing from sub- 
lattice a to the center of the tetrahedron, the hopping parame- 
ters 9 and 9' are transformed to -9 and -9' respectively. In 
other words, the Hamiltonian parametrized by hopping pa- 
rameters (f, 9), (f , 0', 4>') yields identical features as the Hamil- 
tonian parametrized by {t, —9), [f, —ff, cf)'). 

We can also derive hopping matrices obeying the above 
symmetry constraints by considering the Slater-Koster ap- 
proximation of orbital overlaps. In terms of this microscopic 
approach, the relevant Slater-Koster parameters are: (1) Ir- 
Ir hopping via direct overlap of d-orbitals. There are three 
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are of second order. 

At large C/'s, the system is gapped, while for U values 
near the onset of magnetic order, the single-particle spectrum 
may continue to remain gapless even after the onset of mag- 
netic order. The topological Weyl semimetal (TWS) and the 
magnetically-ordered metallic (mAIAO) phases are realized 
in this gapless window — the former is developed via the split- 
ting of the quadratic-band touching at the Fermi-level of the 
SM phase while the latter is realized due to the presence of 
particle-hole pockets in the M phase. 



IV. MAGNETIC EXCITATION SPECTRUM 



FIG. 2. (Color online) Mean-field phase diagram with hopping am- 
plitudes 0' = 0, f'/r = 0.1, f = 1, and f = 57r/6. The phases in 
the non-interacting limit is shown at f/ = 0. First order transitions 
are indicated by dashed lines while second order transitions are in- 
dicated by solid black lines. Lines A (0 = 1.45) and B (6* = 1.70) 
represent the two cuts for which the RPA dynamic structure factors 
are calculated (Fig. 3). 



such overlaps: to- (for cr-bonds), t„ (for 7r-bonds), and ts (for 
5-bonds) for NN (and primed ones for NNN), and (2) to de- 
scribing the hopping between the Ir atoms via the intermedi- 
ate oxygen. Details of such a derivation can be found in Ref. 
3, 6, and 34 and Appendix A. 



III. MEAN-FIELD PHASE DIAGRAM 



As pointed out in Sec. Ill, for a sufficiently large Hub- 
bard interaction U and for 9,12 < 9 < n - 6,12, the mean- 
field solution is the AIAO state. We want to study the na- 
ture of the low-energy magnetic excitation spectrum of this 
magnetically-ordered state, which are composed of the trans- 
verse fluctuations of the spins about their local ordering direc- 
tions. We study these excitations in both the intermediate and 
strong electron correlation regimes and we accomplish this by 
applying two different and contrasting approaches. For the 
case of intermediate- f/, we study the spin- waves by comput- 
ing the RPA transverse spin-spin dynamic structure factor at 
zero temperature. In the large- [/ regime, we perform a strong- 
coupling expansion of our Hubbard model (Eq. 1) to derive 
an effective spin-model. The spin-model is then used to calcu- 
late the spin-wave spectrum about the AIAO state within the 
Holstein-Primakoff approximation. 



In Fig. 2, we show the mean-field phase diagram for 
< 6* < ;r. We have chosen 9 = 9', t'/t = 0.1 and </>' = 57r/6. 
Throughout our calculation we have set f = 1 as our energy 
scale and we comment on its possible values in physical sys- 
tems in the concluding section. We note that this parametriza- 
tion is somewhat different from those used in Ref. 3 where 
the Ir-O-Ir hopping amplitude was set as the unit of energy. 
The phase diagram for > 9 > -n will be identical due to 
the (f, 9), (f', 9', (p') ~ (t, -9), (f', -9', <p') structure as discussed 
in the Section II. Our phase diagram is consistent with those 
obtained in previous calculations. 

In the non-interacting limit, depending on the value of 9, 
we find a strong topological-insulator (STI), the metallic (M), 
and the semimetallic phase (SM). The M phase has small 
particle- and hole-like pockets, whereas the SM phase has a 
quadratic-band-touching at the F-point at the chemical poten- 
tial. At finite U, only two magnetic configurations are found: 
the all-in/all-out (AIAO) and the rotated all-in/all-out (AIAO') 
orders. The AIAO configuration is realized by increasing U 
starting with the SM or M phase. On the other hand, the 
AIAO' configuration is realized by increasing U in the M or 
STI phase. Phase transitions to the AIAO' by an increase in 
U is of first order Also, the transitions between the AIAO 
and the AIAO' phases are of first order and occur at 9,12 and 
n — 9,12. In the non-interacting limit, band-inversion at the 
F-point also occurs at these values of 9. All other transitions 



A. Intermediate-U: RPA dynamic structure factor 

The information about the spin-waves are contained in the 
transverse part of the RPA dynamic spin-spin susceptibility 
matrix. This is given by: 

A'RPA±(q, w) = (1 - UxM¥±{f{, t^))"'OrMF±(q, w)), (9) 

where xmfi. is the bare mean-field transverse spin-spin sus- 
ceptibility, and U is the Hubbard repulsion. Since the py- 
rochlore unit cell has four sublattices, j^^mfj. is a 4 x 4 matrix. 
We compute the trace of the imaginary part of the RPA sus- 
ceptibility, i.e., the RPA dynamic structure factor This trace 
sums over the contribution of the individual spin-wave bands 
(there are four such bands) and gives the overall intensity that 
will be observed in inelastic neutron scattering or RIXS ex- 
periments. The details of the form of the susceptibility matrix 
is discussed in Appendix B. 

Results: We consider the RPA dynamic structure factor 
along two representative cuts (A and B) in the phase diagram 
as U is increased. For cut A,9 - \ .45 while for cut B, = 1 .70 
as shown in Fig. 2. Due to the presence of spin-orbit coupling 
(f2 0), spin-rotation symmetry is explicitly broken, and the 
spin-wave spectra are expected to be gapped. 

The dynamic structure factors for cut A are shown in Fig. 
3a. For this cut, the system is in the quadratic-band-touching 
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(a) (b) 

FIG. 3. (Color online) RPA dynamic structure factors for (a) 9 = 1.45 and (b) = 1.70 for various values of f/; 6oi is in units of f = 1. These 
correspond to cuts A and B in Fig. 2 respectively. Sharp dispersions can be seen at larger values of f/'s, while lighter intensity and broadened 
spectra is seen due to Landau damping in parts of the Brillouin zone at lower U's. Degeneracies at F, X, and W are symmetry-protected and 
are robust features of the AIAO state. 



SM phase in the non-interacting hmit. For U - 4.5, the low- 
lying particle-hole continuum damps and broadens the spin- 
wave spectrum throughout most of the Brillouin zone. This 
so-called Landau damping occurs because spin-waves decay 
through interactions with single-particle excitations and ac- 
quire a finite lifetime, which broadens its spectrum^^^""*^. The 
damping does not occur near the L-point, where the disper- 
sions extend out of the continuum and produce sharp bands. 
As U is increased, the particle-hole continuum is shifted up- 
wards in energy throughout the Brillouin zone, revealing all 
four spin- wave modes. The lower energy modes are relatively 
dispersionless compared to the higher energy modes, which 
disperses most markedly near the L-point. At [/ = 6.5, the de- 
generacies of the spectrum become more apparent: there are 
two two-fold degeneracies at the X- and W-points, while the 
F-point has a three-fold degeneracy. We note that for U < 4.5, 
which includes the TWS phase, all the spin-wave modes are 
damped by the particle-hole continuum. 

For the second cut (Fig. 3b), the system is in the metal- 
lic phase in the non-interacting limit. Well-defined spin- 
wave excitations are only observed for U > 4.0. Starting 
at U - 4.0, the spin-wave modes appear near the F point. 



At U - 4.5, low-energy, damped features can be seen near 
the L-point and along the F - K line. For U = 5.0, most of 
the spin-wave modes become sharply defined as the particle- 
continuum shifts upward. As U is increased further, the low- 
lying dispersion at L shifts up and bands at F begin to sepa- 
rate while maintaining the three-fold degeneracy required by 
symmetry. At {/ = 6.0, the spectrum begins to resemble the 
U - 6.5 spectrum of the first cut. 

Degeneracies at the high symmetry points F, X, and W are 
symmetry-protected, and therefore, they are characteristic to 
the AIAO state, which preserves the lattice symmetry. The 
almost flat dispersion encountered at the zone boundary (X- 
W) is also a distinguishing feature. 



B. Strong-Coupling Expansion - Linear Spin- Wave Theory 

We now look at the spin-wave spectrum in the strong- 
coupling limit of large U/t. In this limit and at half-filling, 
we can apply perturbation theory to obtain the following ef- 
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fective spin Hamiltonian at the lowest order. ^'^ 



^Y.iJSi-Sj + ^u-SixSj + sirpfi, (10) 



where the three terms in the last line are the trace, trace- 
less antisymmetric, and traceless symmetric parts of A)*. 
These terms correspond to the Heisenberg, the Dzyaloshinski- 
Moriya (DM), and the anisotropic interactions, respectively, 
and are related to the hopping amplitudes of Eq. ( 1) by: 

7 = 4C/-'(f2-|v|-/3), 

rff = 8f/-'(v,/v,/-5"W/3), (11) 

(the magnitude of v,/= |v|) is site independent) which holds 
for both NN and NNN hopping amplitudes. The "direct" con- 
figuration of the DM vectors are known to stabilize the AIAO 
state-'^ which is in agreement with our earlier mean-field phase 
diagram. Hence we consider the low-energy spin-wave expan- 
sion about the AIAO state for the above spin Hamiltonian. 

To obtain the spin-wave expansion about the AIAO state 
which orders non-collinearly, we rotate our spin quantization 
axis locally in alignment with the magnetic ordering. To this 
end, we define rotated spin operators, S, such that their local 
5- points to the direction of magnetic ordering at that site. 



5f = (R,(S,))"=/?f5f, 



(12) 



where R,- and R"'' are the rotation operator and its matrix rep- 
resentation that takes the direction of magnetic order at site 
/ and rotates it to the z-axis of the global coordinate system. 
With these rotated operators, we can rewrite Eq. (10): 



(13) 



where R^ indicates matrix transposition. 

After recasting our spin operators in the rotated coordinate 
system, we are in the position to analyze the spin-waves about 
the AIAO state by applying linear spin-wave theory. First, we 
rewrite our spin operators in the Holstein-Primakoff bosonic 
representation: 



(14) 



SJ^ - a] sjls — a\ai, 
S] - s - a'.ai. 



where s is the total spin angular momentum and we have in- 
troduced four flavors of bosons, one for each sublattice of the 
pyrochlore unit cell. Next, we expand and truncate the spin 
Hamiltonian to quadratic order, Fourier transform the bosonic 
operators, and solve for the resulting excitation spectrum via 
a Bogoliubov transformation. 
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FIG. 4. Evolution of the Holstein-Primakoff spin-wave spectra as 
angle 6 is varied between (e,/2 « 0.96) < 9 < (n - 9,/2 « 2.19); oj 
is in units of f/U. 6 and n - 9 yield identical spin- wave spectra and 
this redundancy is reflected by each plot having two 9 values. Like 
the RPA results in Fig. 3, degeneracies at F, X, and W are symmetry- 
protected. 



Results: We consider the spin-wave spectra obtained from 
the effective spin Hamiltonian generated by only NN hopping 
amplitudes. Adding up to f' - 0.1 f NNN hopping ampli- 
tudes (not shown) only leads to small changes (see below). 
In Fig. 4, we depict the evolution of the spin-wave spectrum 
as the NN hopping parameter 6 is varied. In the absence of 
NNN hopping amplitudes, 6 and n - yield identical spin- 
wave spectra. This structure is again related to our choice 
of angular parametrization: and -9 are related by a basis 
transformation and should therefore have the same spin-wave 
spectrum. Moreover, -0 ^ -0 + nis equivalent to r — > -f and 
V -V, which leaves J, D,j, and invariant. Hence, and 
7T-0 yield the same spin-wave spectrum and this fact is noted 
by assigning each of the plots in Fig. 4 with two values. 

For 0i/2 < < TT - 0,/2, the spin-wave spectrum is gapped 
and the gap decreases as we approach either endpoints of the 
interval. At the endpoints, two of the four bands of the spec- 
trum become both gapless and dispersionless, while outside 
the endpoints, the lowest bands become negative in energy, 
signaling an instability of the AIAO state. The onset of these 
instabilities is consistent with NN mean-field theory results, 
which predicts first order transitions (between the AIAO to 
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the AIAO' phase) at = 6,12 and = tt - 6,12. 

We note that the degeneracies at the high-symmetry points 
r, X, and W are consistent with the RPA results as they are 
protected by symmetry. Also, the flat dispersions at the zone 
boundary (X-W) is also encountered in the present spin-wave 
calculation. On the other hand, the lowest energy dispersion 
along the L-F line is absolutely flat in this NN model. How- 
ever, on adding small NNN hopping (up to f =0.1 t), it ac- 
quires small dispersion. This should be contrasted with RPA 
results in Fig. 3 where modes along the L-F line are more 
strongly dispersive. 



C. Comparison of RPA and strong-coupling results 

As the ratio of typical hopping scale to the Hubbard re- 
pulsion scale (tjU) increases, the higher-order contributions 
to the strong-coupling perturbative expansion (in tjU) be- 
come increasingly important and the strictly NN model we 
employed in Section IV B becomes inadequate in describ- 
ing the magnetic excitations of the AIAO state. Not only do 
higher-order contributions generate further-neighbor Heisen- 
berg exchanges, ring-exchange type terms arise and lead to 
renormalization of NN quadratic terms at the linear spin-wave 
level. Therefore, we should not expect a perfect agreement 
between the RPA and strong-coupling results. 

Nevertheless, we attempt to fit (by eye estimation) the RPA 
results for large U {e.g., U — 6.0 in cut B) with a linear spin- 
wave (LSW) spectrum. First, we fit the RPA dispersion fea- 
tures and overall bandwidth, resulting in 7 = 0.13, |D| =0.13, 
|F| = 0.07 (see Eq. 1 1 for definitions). These parameter val- 
ues are diff'erent from those obtained from Eq. 1 1 which is 
based strictly on the leading order strong-coupling expansion 
for the NNs. The gap obtained from the RPA calculation is 
Arpa ^ 1.70. 

The resulting LSW fit captures the dispersion along F-X- 
W quite well, but fails to capture the low-lying modes along 
the L-F line and, in general, the higher energy modes where 
the fit at best is qualitative. We would also like to point out 
that fits at lower U and along cut A have been attempted but 
large discrepancies in both the features of the spectrum and 
the spin-wave gap have been found. 

The above fitting results points out that the NN spin-model 
(and also the NNN spin-model with up to 10% NNN hop- 
ping amplitude) is grossly inadequate to fit the RPA spin- 
wave spectrum quantitatively for parameter values that en- 
compasses the regime appropriate for the pyrochlore iridates. 
In this context we would like to point out that recent estimates 
of small charge gap (~ 10 meV) from the resistivity measure- 
ments in EuiIraOv seems to suggest that the intermediate- 
coupling calculations may be better suited to describe this 
compound. 



V. SUMMARY 

To summarize, we have calculated the structure of the mag- 
netic excitation spectrum for the AIAO state that has been pro- 
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FIG. 5. (Color online) RPA results at t/ = 6.0 for cut B (blue and 
yellow color-map) overlaid with a fit with the NN linear spin-wave 
spectrum (red); oj is in units of / = 1. Though some features along 
the F-X-W line can be fitted, other high-symmetry lines show a larger 
discrepancy. 



posed for the pyrochlore iridates such as Eu2lr207 or Yt^t^iO-j. 
For intermediate correlations, we have calculated the trans- 
verse spin-spin dynamic structure factor within the RPA ap- 
proximation. Features particular to the AIAO configuration 
that can lead to conclusive identification of the magnetic or- 
der in Eu2lr207 and Yii^iOi were discussed. For the large- [/ 
limit, we used a strong-coupling expansion to derive a spin- 
model and calculated the linear spin-wave spectrum using the 
Holstein-Primakoff approximation. By fitting the RPA results 
with the linear spin-wave spectrum, we showed that results in 
the intermediate correlation regime are substantially different 
from those obtained from the strong-coupling theory. 

Finally, from our calculations, we can make rough esti- 
mates of the experimental energy scales for the excitation gap 
and dispersion bandwidth by using Slater- Koster parametriza- 
tion of orbital overlaps (the relation connecting the symmetry- 
allowed parameters in Eq. 6 to Slater-Koster parameters are 
detailed in Appendix A). Similar to Refs. 3 and 34, we have 
used t„ = -2ta-/3, ts = 0, while = -1.16fo(-1.37fo) 
for cut A{E) (the different orbital overlaps were introduced 
towards the end of Section II and are also used in Ap- 
pendix A). Typically, the value of t,, — the Ir-O-Ir hopping — 
is about 200-350 meV in iridates with octahedral oxygen 
environments.'*"'*' Here we choose a representative value of 
to ~ 300 meV. With these values, in Fig. 3, the spin gap found 
from our RPA calculation is on the order of 100 meV and the 
dispersion width is on the order of ~ 15 meV. This value of 
spin gap is found to be very sensitive to the value of U while 
the bandwidth always remains in the same regime. While a 
100 meV gap can be resolved within current RIXS resolution, 
the dispersion width may be presently on the borderline of re- 
solvability. 

In other iridium compounds, recent progress have been 
made in RIXS^^"'-"'"*", neutron scattering^^-''*^', resonant mag- 
netic x-ray scattering (RMXS)"*"* '*'', and x-ray absorption spec- 
troscopy (XAS)'*^ experiments. In particular, RIXS and in- 
elastic neutron scattering experiments have recently revealed 
the magnetic excitation spectra of Sr2lr04^", Sr3lr207^^', and 
Na2lr03^-. Future applications of these techniques and im- 
provements in experimental resolution may help reveal the 
magnetic behavior and, in particular, the magnetic excitation 
spectra of pyrochlore iridates, thereby conclusively determin- 
ing the nature of their magnetic order. 
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parameters is given by the following relations 



Appendix A: Microscopic origins of hopping parameters 



V2 



1 



8748 

V2 

8748 
1 

8748 



316f;,-43f<s + 520fo), 
l6Qt„ - 22Qts + lUta) , 
(6994- 16284- 1843?^), 

(-1564-27204-4?^), 
(-1444-9604+ 1104f9, 



(A2) 



From these relations, it is straightforward to use Eqs. 7 and 
8 to relate the angular parameters to the above Slater-Koster 
parameters. 

The range of physical NN hopping parameters explored in 
Ref. 34 (-1.2 < < -0.5, = -2f^/3, = and = 1) 
corresponds to 0.85 < < 1.51 with the appropriate energy 
scaling of t. 



Appendix B: Mean-field transverse spin-spin susceptibility 



In order to estimate the values of the hopping parameters, 
we turn to a microscopic analysis of hopping paths. We briefly 
discuss the results of such an analysis and refer to Ref. 5, 6, 
and 34 for more details. 

We consider two types of hopping: Ir-Ir hopping via over- 
lap of (i-orbitals and O-Ir hopping between p-orbitals of O 
and c/-orbitals of Ir We parametrize d-d overlaps with Slater- 
Koster amplitudes tg-, t„, and tg (and primed ones for NNN), 
while for p-d overlaps, we parametrize them with amplitudes 
tpdcr, tpdn, and O-Ir occupation energy difference e. Here the 
subscripts cr, tt and S denotes the type of overlap of the or- 
bitals. Also, in this microscopic picture, we always work in 
the local axes defined by the oxygen octahedra surrounding 
each Ir 

To arrive at a 7efF - 1/2 model, we first employ second 
order perturbation on the O-Ir hopping to generate an effective 
NN Ir-Ir hopping between c/-orbitals. This indirect hopping is 
given by 



= tj\e\ (Al) 



where e is the diff'erence of the on-site charging energies be- 
tween the oxygen 2p orbitals and Ir 5d orbitals. (This indi- 
rect hopping receives contribution from the tpda- overlap in 
presence of distortion in the oxygen octahedra). We now 
project the t/-orbitals into the local t2g and, finally, into the 
local 7eff - 1/2 basis, to find the hopping matrix which has 
the form given by Eq. (6), where the relation with the effec- 
tive hopping parameters and more microscopic Slater-Koster 



As Eq. 9 indicates, we need to compute the mean-field 
transverse spin-spin susceptibility in order to obtain the RPA 
susceptibility. Here we provide details of such calculation for 
non-collinear magnetic order like the AIAO state. 

The mean-field susceptibility matrix is given by: 

Xmf±(^' = -'0(0 ([(Sl(q, 0)/, (Sl(-q, 0)),]) , (Bl) 

where S"(q, f) denotes the Fourier transform of the compo- 
nent of the spin operator that is perpendicular to the magnetic 
ordering direction for sublattice a (see Fig. 1 for sublattice 
convention), and ; indexes the components of S" (q, f), which 
are to be summed over. If the direction of magnetic moment 
on sublattice is given by the unit vector n", then the transverse 
spin operator for that sublattice is given by 

S^(q) = -h" X (n" X S"(q)), (B2) 

where S''(q) is the Fourier transform of the spin operator. For 
the AIAO state, h" points along the local [111] direction of the 
lattice. 

We write the spin operator with electron operators as 

S"(q)=^2<.q,«'^''/^^M (B3) 
^ k 

where a,/3 -1, i. Using this we get 
X 2 {^"((1 + ki, f)Cfl/j(^i, t), cl^(-q + k2, 0)cM(k2, 0)]) , 

k,k2 

with 

KyS = 5 («° X («" X ^afi)) ■ 5 («' X [h" X CT,,)) , (B5) 

Lastly, we transform our basis to the band basis using the 
results from our Hartree-Fock mean-field calculation, evalu- 
ate the two-body expectation value via Wick's theorem, and 
Fourier transform to frequency space to obtain ;^f^p^(q, a>). 
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